Coordination motifs and large-scale structural organization in atomic clusters 
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The structure of nanoclusters is complex to describe due to their noncrystalhnity, even though 
bonding and packing constraints limit the local atomic arrangements to only a few types. A compu- 
tational scheme is presented to extract coordination motifs from sample atomic configurations. The 
method is based on a clustering analysis of multipole moments for atoms in the first coodination 
shell. Its power to capture large-scale structural properties is demonstrated by scanning through the 
ground state of the Lennard- Jones and Ceo clusters collected at the Cambridge Cluster Database. 
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I. INTRODUCTION 

Nanoscale atomic clusters have attracted a great deal 
of attention in recent years due to their promising ap- 
plications in catalysis, photonics, and bioimagingi^i^ 
The functionality of these clusters is controlled by their 
structure, which may assume a variety of forms as a re- 
sult of the generally complex energy landscape governing 
the low energy cluster conformations."* Atomic arrange- 
ments in these clusters generally do not follow the crystal 
structure of the corresponding bulk material, but may 
nevertheless possess certain regularities. Suitable repre- 
sentation of these structures, particularly for small and 
medium sized clusters, is important for uncovering mech- 
anisms that underly their thermodynamic stability and 
functionality. 

Traditionally, atomic clusters of high point symme- 
try are described in terms of the five Platonic solids 
and their truncated polyhedra.-i^ Starting from, say, a 
13-atom core as illustrated in Fig. 1, successive layers 
of atoms can be added to produce clusters of "magic 
size" which feature prominently in spectroscopic mea- 
surements. Two particularly well-known series are the 
Mackay icosahedra^ (MIC) that grow from Fig. 1(c) 
and the Marks decahedral (MD) that grow from Fig. 
1(d). Noble-gas nanoclusters follow the Mackay icosahe- 
dra series up to a few thousand atoms before switching 
to the Marks decahedra series.'* Both types of clusters 
can be viewed as multiply twined particles that take ad- 
vantage of the low energy surfaces of the constituting 
face-centered-cubic (fee) grains. At even larger sizes, ac- 
cumulation of strain in the particle drives a transition to 
the single crystal structure, although the size at which 
such transition takes place is still debatedi^ 

Real structures often contain defects and imperfections 
as compared to the ideal ones Recently, Polak et al. in- 
vestigated medium-sized Lennard- Jones (LJ) clusters of 
a few hundred atoms obtained in Monte Carlo simula- 
tions at low temperatures f^ii^ They proposed a scheme 
to describe the interior structure of an atomic cluster in 
terms of the four "coordination polyhedra" of Fig. 1. In- 
deed, when atoms with a particular type of coordination 
are displayed, large-scale features such as cubic domains 
or polyicosahedral grains become readily visible. Polak 



FIG. 1: (Color online) The first coordination shell at co- 
ordination number 12. (a) face-centered-cubic (fee); (b) 
hexagonal-close-packed (hep); (c) icosahedral (ico); (d) deca- 
hedral (dec). 



applied the scheme to describe the structure evolution 
during cooling runs, and in cluster growth starting from 
a seed crystallite. It is conceivable that insights gener- 
ated from such studies can assist experimental designs 
aimed at producing a particular class of nanoclusters. 

The four coordination polyhedra in Fig. 1 were iden- 
tified by Polak through visual inspection of a large num- 
ber of atomic configurations. Here we extend Polak's 
work to construct a systematic and robust computational 
methodology for local coordination analysis and large- 
scale structure identification. This is necessary as the 
local atomic structures often contain distortions or even 
belong to unknown categories. The idea is quite sim- 
ple. Given a set of atomic coordinates, we first compute 
the radial distribution function to determine a suitable 
cutoff distance for the first coordination shell. We then 
parametrize the coordination shell configuration of each 
interior atom using the multipole moments Qi of the 
shell atoms." The data set is examined for their clus- 
tering properties in the "shape space" spanned by the 
Qis. Each well-defined cluster in the shape space rep- 
resents a re-occuring local structural unit which we call 
a "coordination motif". Atoms in the same cluster are 
then assigned the same "type". The scheme allows for 
fast classification of interior atoms and their local neigh- 
borhoods. The geometric shape of individual motifs and 
their spatial organization provide valuable information 
on large-scale structure organization. The MIC and the 
MD series, for example, each possesses a unique composi- 
tion profile of the coordination motifs. In such cases, the 
motif content alone can be used for accurate structure 
association. 
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The paper is organized as follows. In Sec. II we de- 
scribe the computational procedure for the identification 
of structure motifs. The ground state LJ clusters are 
used as examples. Section III contains a discussion of 
large-scale structure organization based on the spatial 
distribution of motifs of different types. Particular em- 
phasis is given to the decagonal motifs that define local 
pentagonal symmetry axis for twinning. Section VI con- 
tains a summary of our findings and conclusions from the 
work. 



II. LOCAL STRUCTURE MOTIF 
IDENTIFICATION 

The procedure for the coordination motif identification 
consists of three steps: i) specification of the first coordi- 
nation shell through an analysis of the radial distribution 
function; ii) parametrization of the atomic configuration 
for completed shells; iii) clustering of the data to deter- 
mine the most frequently occuring coordination motifs. 

We illustrate the method in the description of the (pu- 
tative) ground state of LJ clusters with up to N = 1610 
atoms, available at the Cambridge Cluster Database 
(CCD).^^'^^ Figure 2(a) shows the radial distribution 
function obtained from the histogram of pairwise dis- 
tances rij = \ri — Vjl for three selected atomic clus- 
ters at iV = 55, 549 (complete Mackay icosahedra) and 
1200, at a bin size Ar — 0.01 (in units of the LJ ra- 
dius a). Note that the minimum of the LJ potential is 
at Tniin = 2^/^CT ~ 1.12a. In all three cases, there are 
well-defined distance gaps which allow for unambiguous 
identification of first and second coordination shells for 
all atoms on the cluster. The position of the first two 
peaks are at ri ~ 1.1 and r2 — 1.55, respectively, with 
the ratio ^2/^1 ~ 1.4 which corresponds well to the value 
1/2 for the fee or hep structure. The plot suggests a cut- 
off distance Tcut = 1-3 for atoms to be included in the 
first coordination shell. To distinguish interior from sur- 
face atoms, we compute the offset distance roSset between 
the center atom and the center of mass (CM) of its first 
coordination shell. Figure 2(b) shows a scatter plot of 
''offset against the coordination number for all atoms on 
the TV = 549 cluster. The approximately linear behav- 
ior of the data indicates that completed shells have 12 
atoms. 

The 12-atom neighborhoods as shown in Fig. 1 are 
distinguished from one another by the angular distribu- 
tion of atoms on the shell. Provided the distance of shell 
atoms to the center atom falls in a narrow range, the 
shell configuration is well represented by a density func- 
tion p{6^(j)) = X]fc^(^ ~ (^kT(t> ~ 4>k) on the unit sphere, 
where {0k,4'k) are the polar coordinates of shell atom k 
with the origin on the center atom. Consider now the 
"multipole" expansion in the spherical harmonics. 
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FIG. 2: (Color online) (a) Radial distribution function of LJ 
clusters at three selected sizes, (b) The distance between the 
center of mass (CM) of the first coordination shell and the 
center atom against the coodination number for all atoms in 
a. N = 549 cluster. Here rcut = 1.3. 



where a/,„ = J2k ^imi^k.^ 4>k)- The multipole moments of 
the shell configuration are given hy^ 
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where C is the number of shell atoms, and the over line 
bar denotes average over the shell. Note that Qi = l (all 
X) if there is only one atom on the shell. It is easy to 
verify that the multipole moments Qi are invariant with 
respect to a rigid rotation of the shell density. 

In Ref.^^, the moments and Qg are used to dis- 
tinguish the four atomic configurations shown in Fig. 1. 
We propose here a more general approach based on the 
clustering property of interior atoms in the shape space 
spanned by the multipole moments \Qi\. Atoms that are 
close to each other in this representation have very sim- 
ilar neighborhoods. Strong clustering of the data points 
in the shape space suggests existence of only a few coor- 
dination motifs. Standard clustering methods^ can then 
be used to assign coordination motifs to the atoms. On 
the other hand, weak clustering would imply a more con- 
tinuous spectrum of local neighborhoods without a clear 
local structural preference. 

Figure 3 shows the multipole moments (54, and Qg 
for all interior atoms on the ground-state LJ clusters at 
the three selected sizes. In all cases, the data exhibit 
very strong clustering which agrees well to the four shell 
configurations in Fig. 1. Note that the cluster that corre- 
sponds to the decagonal motif has a slightly bigger spread 
than the other three clusters. Inspection of the atomic 
coordinates show that the decahedron has a shorter axis 
in the two smaller systems and is elongated slightly in 
the iV = 1200 case. The moment Qg performs better 
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FIG. 3: (Color online) Scatter plot of the multipole moments 
of interior atoms on the ground state LJ clusters at three 
different cluster sizes, (a) Qe versus Q4; (b) Qs versus Qs- 
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than Q4 in discriminating the hep and dec motifs. 

Thermal fluctuations lead to distortions of the polyhe- 
dra so that the corresponding data points in the shape 
space are more spread out than in the ground state. We 
have investigated the dependence of the width of such 
clouds on the noise strength using simulated data. Start- 
ing from each of the four regular polyhedra, we randomly 
displace the 12 shell atoms laterally to create noisy den- 
sity patterns. The displacement of each atom is restricted 
to within a distance Ar = RA9 from its original position, 
were R is the radius of the coordination shell and Ad de- 
fines the noise strength. Figures 4(a)- (c) show a compar- 
ison of the Qi's for the ideal (sohd symbols) and noisy 
(open symbols) structures at A6 — O.OTtt. In the latter 
case, 100 shell configurations are generated from each of 
the four regular polyhedra. Error bars on the data in Fig. 
4(b) represent standard deviations over the 100 samples. 
From the plots it is seen that the icosahedral coordination 
motif is well separated from the other three and can be 
easily recognized. The hep motif is sandwiched between 
the fee and dec motifs in the shape space and is more 
easily deformable into either of the two. At A9 = O.OTtt, 
the hep cluster already begin to merge with the fee and 
dec clusters. The data (with error bars) in Fig. 4(b) 
suggest that either Qq or can be used to separate the 
hep from the fee motif while Qg can be used to separate 
the hep from the dec motif. Enhanced performance is 
expected when the classification is done in the high di- 
mensional shape space which combines information from 
all the Qi's. 

We have also generated noisy atomic clusters by adding 
random displacements to the ground state LJ configura- 
tions. Each component of the atomic displacement is 
taken to be gaussian distributed with zero mean and 
standard deviation h. The resulting atomic cluster is an- 
alyzed following the procedure described above. Figure 
4(d) shows an example of the local coordination data ob- 
tained from 10 randomized samples of the N = 549 clus- 



FIG. 4: (Color online) The multipole moments Qi for regular 
(solid symbols) and noisy (open symbols) polyhedra. (a) Qi 
against I in the noiseless case; (b) Qi against I when shell 
atoms are displaced laterally at a noise strength A6 — O.OTtt; 
(c) Qe versus Qg for individual shell configurations at A6 = 
O.OTtt; (d) Qe versus Qg from 10 randomized samples of the 
N = 549 LJ cluster at Ar = O.lrmin (see text). 

ter at h = O.OGrmin or Ar = S^^^h = O.lrmin- Only atoms 
with coordination number 12 are included. The symbols 
used for the selected atoms in the plot correspond to their 
coordination type in the ground state cluster. It is seen 
that the data at this noise level still form well-defined 
clusters in the shape space, with the majority of atoms 
keep their original identity. There are however a few out- 
liers (points with Qg less than 0.35) where the shell con- 
figuration has changed so significantly that they can no 
longer be classified into any of the four categories. Closer 
inspection of the atomic configurations indicate that such 
cases arc usually accompanied by large displacement of 
atoms in the neighborhood so that a more careful defini- 
tion of the first coordination shell for the atom in ques- 
tion is needed. Interestingly, the Lindemann criterion 
for crystal melting sets a limit on Ar to be around 10- 
12% of the nearest neighbor distance at the bulk melting 
point.~ From our simulation studies, this is also the noise 
level around which different coordination motifs start to 
become indistinguishable from each other. 



III. COMPOSITION PROFILE AND 
LARGE-SCALE STRUCTURE ORGANIZATION 
OF ATOMIC CLUSTERS 

Using the moment pair (Qg, Qg) as classifier, we com- 
puted the composition profile Cq — Na/NintcrioT of the 
four coordination motifs for all ground state LJ clusters 



TABLE I: Number of icosahedral (ico), decahedral (dec), 
hexagonal-close-packed (hep) and face-centered-cubic (fee) 
coordinated atoms in the complete Mackay icosahedra (MIC) 
and Marks decahedra (MD) series. The total number of com- 
pleted shells (including surface atoms) is given by n. 
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in the size range = 13 - 1610. Here Na is the number 
of a = fee, hep, ico, dec atoms in the cluster interior 
and iVintorior is the total number of interior atoms. Re- 
sults are shown in Fig. 5. Each curve gives the percent- 
age of a particular motif in the atomic cluster against the 
cluster size. As a reference, Table I lists the number of 
interior atoms of each type for complete Mackay icosahe- 
dra and Marks decahedra as a function of the layer num- 
ber n. Also given are the total number of interior atoms 
Mntorior as well as the cluster size N. Using these expres- 
sions, we computed the concentration of the fee motif for 
the two series, as indicated by the dashed and dotted 
lines in Fig. 5. The overall trend of structure evolu- 
tion is evident. Beyond TV = 147 (three completed shells 
in an icosahedral structure), the ground state essentially 
alternates between the two structures. The main transi- 
tion between the icosahedra-dominated regime at smaller 
sizes to the decahedra-dominated regime at larger sizes 
takes place at around N = 1030, indicated by an up- 
ward jump in the Cfcc curve (black solid line) on the plot. 
Windows exist on either side of the transition where the 
minority series pops in. The physical mechanism behind 
these transitions, which has to do with the competition 
between strain and surface energies, has already been 
discussed extensively in the literature^ 

We have examined the spatial distribution of the coor- 
dination motifs on large atomic clusters. Figure 6 shows 
the position of ico (blue), dec (green), and hep (orange) 
atoms at TV = 923 (MIC) and 1103 (MD). The icosahe- 
drally coordinated atom is only present at the center of 
the MIC cluster. The decagonally coordinated atoms are 
located on the five-fold symmetry axes in both clusters, 
while the hep atoms form two-dimensional sheets that 
separate the fee grains (not shown) at symmetry related 
crystallographic orientations. 

Clusters of size less than three atomic layers exhibit 
a greater variety of organizational behavior due to the 
more subtle and delicate balance of packing, strain, and 
surface effects. Figure 7 shows the composition profiles 
in this size range for two representative classes of ground 
states in the CCD: the LJ clusters discussed above and 
the Ceo clusters computed using the Pacheco and Prates- 




FIG. 5: (Color online) Composition profile of four coordi- 
nation motifs in the ground state LJ clusters. Dashed and 
dotted lines give the theoretical values of Cfcc for complete 
Mackay icosahedra and Marks decahedra, respectively. 




FIG. 6: (Color online) Position of ico (blue), dec (green), and 
hep (orange) atoms in the ground state LJ cluster at size (a) 
iV = 923 (Mackay icosahedron) and (b) iV = 1103 (Marks 
decahedron), showing large-scale structural organization. 



Ramalho potentialJ^ Most of the LJ clusters in this size 
range has a single icosahedron motif (open blue circle) at 
the core. Their structure essentially follows that of the 
Mackay icosahedra series, with modifications close to the 
surface layer. Notable exceptions are found at TV = 38 
(fee) and TV = 75 — 77 (Marks decahedron), where the 
icosahedron core is missing. 

The Cgo clusters behave quite differently from the LJ 
clusters. The icosahedron core is present only at TV < 15. 
About one third of the clusters in the size range 31 — 100 
are cubic (i.e., with only fee and hep motifs), whereas 
the rest contain some form of the decahedral structure. 
It has been suggested that, due to the narrow range of 
the molecular potential, Ceo clusters prefer close-packed 
structures as compared to the LJ clusters J^^^S, This is in 
good agreement with the much higher percentage of fee, 
hep, and dec motifs in the composition profiles as shown 
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FIG. 7: (Color online) Motif composition profile for Lennard- 
Jones (LJ) and Ceo clusters below A'' — 100. 



in Fig. 7. We have also examined the different series of 
the ground state Morse clusters in the CCD that corre- 
spond to a broad range of softness of the pair potentiali^ 
Indeed, as the range of the interaction shortens, a change 
from the icosahedral to decagonal and cubic structures is 
observed. 

The decahedral motif plays a key role in the large-scale 
structural organization due to its unique pentagonal sym- 
metry axis. In both the Mackay icosahedra and Marks 
decahedra series, decahedra motifs stack together to form 
linear chains to propagate the (local) pentagonal sym- 
metry all the way to the cluster surface. Two or more 
pentagonal axes can emanate from an icosahedral mo- 
tif or an incomplete icosahedron at the surface. Chains 
of decahedra have also been observed in less ordered LJ 
clusters by Polak Jii^ The network of decahedral motifs 
(plus their orientation), to a great degree, defines the 
overall structure of these low energy atomic clusters. 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have presented a general method to 
identify re-occuring local structural units in an atomic 
cluster without prior information. The structural units, 
taken here to be atoms on the first coordination shell, are 
represented by a set of multipole moments whose cluster- 
ing properties form the basis for robust and efficient clas- 
sification schemes. Applying the scheme to the ground 
state of the LJ and Cgo clusters, we reconfirmed that the 
fee, hep, ico and dec motifs provide a complete set of 
the 12-atom neighborhood for all interior atoms. For the 
LJ clusters, the motif compositions show a systematic 
dependence on cluster size, in agreement with the ex- 
isting understanding on the structure evolution of these 
clusters. Large-scale structural organization are readily 
revealed when interior atoms are displayed according to 
their coordination type. Under this scheme, the level of 
complexity in structure representation is significantly re- 
duced. This "microstructure" characterization can play 
an important role in the development of coarse-grained 
models that combine geometry and energetics to eluci- 
date the equilibrium and kinetic properties of medium- 
sized atomic clusters at low temperatures. 

The study of noisy atomic clusters allows one to relate 
the size and shape of the cloud surrounding each ideal 
shell structure in the shape space to the actual displace- 
ment of atoms in the physical space. Interestingly, the 
noise level at which the clouds for different coordination 
motifs overlap coincides well with the one given by the 
Lindemann criterion for crystal melting. This observa- 
tion offers a potentially fruitful way to analyze melting 
in super-heated solids^^ in terms of the local coordination 
patterns and their spatial distribution. 

We wish to thank Dr. X.-G. Shao for providing the 
LJ cluster coordinates in the size range 1001-1610. This 
work was supported by the Research Grants Council of 
the HKSAR through grant HKBU 2020/04P. 
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